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We discuss the theoretical and experimental status of the CP violating ratio e'/e. We revise our 
1997 standard-model estimate — based on hadronic matrix elements computed in the chiral quark 
model up to 0(p 4 ) in the chiral expansion — by including an improved statistical analysis of the 
uncertainties and updated determination of the Cabibbo-Kobayashi-Maskawa elements and other 
short-distance parameters. Using normal distributions for the experimental input data we find 

Re e'/e = (2.2 ± 0.8) x 10~ 3 , 

whereas a flat scanning gives 

0.9 x 10" 3 < Re e'/e < 4.8 x 10" 3 . 

Both results are in agreement with the current experimental data. The key element in our estimate is, 
as before, the fit of the A7 =1/2 rule, which allows us to absorb most of the theoretical uncertainties 
in the determination of the model-dependent parameters in the hadronic matrix elements. Our semi- 
phenomenological approach leads to numerical stability against variations of the renormalization 
scale and scheme dependence of the short- and long-distance components. The same dynamical 
mechanism at work in the selection rule also explains the larger value obtained for e'/e with respect 
to other estimates. A coherent picture of K — > it-k decays is thus provided. 



I. INTRODUCTION 

The violation of CP symmetry in the kaon system (for two recent textbooks on the subject see ref. 0) is parame- 
terized in terms of the ratios 
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Eqs. (1.1) can be written as 
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where lo = A2/A0 is the ratio between the isospin 1 — 2 and components of the K — > ttt: amplitudes, the anomalous 
smallness of which is known as the AI =1/2 selection rule ||- The complex parameters e and s' are introduced to 
quantify, respectively, indirect (via K^-Kg mixing) and direct (in the Kj, and Kg decays) CP violation. They are 
measurable quantities, and e has been known to be non vanishing since 1964 0. 
The AS* — 1 effective lagrangian Cw is given by 

C w = -Y.Ciin) Qi(jjL) , (1.3) 
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C i ( t ,) = ^v ud v: 
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In ( |l.4| ), Gi? is the Fermi coupling, the functions and yi(n) are the Wilson coefficients and the Cabibbo- 

Kobayashi-Maskawa (CKM) matrix elements; r = —VtdV t * s /V u dV* s . According to the standard parameterization of 
the CKM matrix, in order to determine e'/e , we only need to consider the yi(fj-) components, which control the 
CP- violating part of the lagrangian. The coefficients yi(fJ-), and Zi(/j,) contain all the dependence of short-distance 
physics, and depend on the t, W, b, c masses, the intrinsic QCD scale Aqcd, the 75-scheme used in the regularization 
and the renormal izat ion scale /1. 

The Qi in eq. (1.3) are the effective four-quark operators obtained in the standard model by integrating out the 
vector bosons and the heavy quarks t, b and c. A convenient and by now standard basis includes the following ten 
operators: 
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where a, (3 denote color indices (a,/3 = 1, . . . ,N C ) and e q are the quark charges (e u = 2/3, id — e s — —1/3). Color 
indices for the color singlet operators are omitted. The labels (V ± A) refer to the Dirac structure 7 M (1 ± 75). 

The various operators originate from different diagrams of the fundamental theory At the tree level, we only have 
the current-current operator Q2 induced by IV-exchange. Switching on QCD, the gluonic correction to tree-level W- 
cxchangc induces Q\. Furthermore, QCD induces the gluon penguin operators Qss- The penguin diagrams induce 
different operators because of the splitting of the color quark (vector) current into a right- and a left-handed part 
and the presence of color octet and singlet currents. Electroweak penguin diagrams — where the exchanged gluon is 
replaced by a photon or a Z-boson and box-like diagrams — induce and also a part of Q 3 . The operators Qs.io 
are induced by the QCD renormaliz atio n of the electroweak penguin operators Q7,9- 

Even though the operators in eq. (1.5) are not all independent, this basis is of particular interest for any numerical 
analysis because it has been extensively used for the the calculation of the Wilson coefficients to the next-to-leading 
order (NLO) order in a s and a e 0, in different renormalization schemes. 

In the standard model, e' can be in principle different from zero because the 3x3 CKM matrix Vij , which appears 
in the weak charged currents of the quark mass eigenstate, is in general complex. On the other hand, in other models 
like the superweak theory the only source of CP violation resides in the K°-K° mixing, and s' vanishes. It is 
therefore of great importance for the discussion of the theoretical implications within the standard model and beyond 
to establish the experimental evidence and precise value of e' . 

The ratio e'/e (for a review see, e.g., p|-||]) is computed as 
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where, the CKM combination X t = VtdVts an d, referring to the AS* = 1 quark lagrangian of eq. (1.3) 
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and (Q i ) I = (2ir,I\Q i \K). 

We take the phase <fi = ir/2 + Sq — 82 — f = (0 ± 4)° as vanishing and we assume everywhere that CPT is 
conserved. Therefore, Re e'/e — e'/e. 



Notice the explicit presence of the final-state- interaction (FSI) phases 5i in eqs. (T7) and (b_8). Their presence is 
a consequen ce o f writing the absolute values of the amplitudes in term of their dispersive parts. Theoretically, given 
that in eq. (1.4) r <C 1, we obtain 



tan(5 7 _ 

Finally, f2, ; + r; ' is the isospin breaking (for m u 7^ mj) contribution of the mixing of 7r with 77 and 77'. 



(1.9) 
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A. Preliminary remarks 



The experiments in the early 1990's [[10 11 could not establish the existence of direct CP violation because they 
agreed only marginally (one of them being consistent with zero [ p"l[ ) and did not have the required accuracy. Dur ing 
1999, as we shall briefly recall in the next section, the preliminary analysis of the new run of experiments (lj H 
have settled on a range of values consistent with the previous NA31 result, conclusively excluding a vanishing e' . On 
the other hand, in order to asses precisely the value we must wait for the completion of the data analysis which will 
improve the present accuracy by a factor 2-3. 

On the theoretical side, progress has been slow as well because of the intrinsic difficulty of a computation that 
spans en ergy scales as different as the pion and the top quark masses. Nevertheless, the estimates available before 
1999 [|l5|-jl7[ pointed to a non- vanishing and positive value, with one of them fl7j| being in the ball park of the present 
experimental result. 

We revise our 1997 estimate 0] of e'/e by updating the values of the short-distance input parameters — among 
which the improved determination of the relevant CKM entries — and by including the Gaussian sampling of the 
experimental input data. We also update the values and ranges of the "long-distance" model parameters (quark and 
gluon condensates, and constituent quark mass), by including a larger theoretical "systematic" error (±30%) in the 
fit of the CP conserving K — > irir amplitudes in order to better account for the error related to the truncation of the 
chiral expansion. 

For the sake of comparison with other approaches, we give our results in terms of the so-called -^-parameters, in 
two 75 regularization schemes: 't Hooft-Veltman (HV) and Naive Dimensional Regularization (NDR). 

The combined effect of the new ranges of the input parameters and ImAt makes the central value of e'/e slightly 
higher than before, while the statistical analysis of the input parameters reduces the final uncertainty. For a more 
conservative assessment of the error, we give the full range of uncertainty obtained by the flat span of the allowed 
ranges of the input parameters. The result is numerically stable as we vary the renormalization scale and scheme. 
We conclude by briefly reviewing other estimates of e' published in the last year. 



II. THE EXPERIMENTAL STATUS 



Experimentally the ratio e'/e is extracted, by collecting and Kg decays into pairs of ir° and tt , from the 
relation 
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and the determination of the ratios rj^ and 7700 given in (1.1). 

The announcement last year of the preliminary result from the KTeV collaboration (FNAL) (l2) 

Re e/e' = (2.8 ± 0.41) x 10~ 3 , 

based on data collected in 1996-97, and the present result from the NA48 collaboration (CERN), 

Re e/e' = (1.4 ±0.43) x 10~ 3 , 



(2.1) 



(2.2) 



(2.3) 



based on data collected in 1997 |l3| and 1998 ]14|] , settle the long-standing issue of the presence of direct CP violation 
in kaon decays. However, a clearcut determination of the actual value of e' at the precision of a few parts in 10 must 
wait for further statistics and scrutiny of the experimental systematics. 

By computing the average among the two 1992 experiments (NA31 and E731 Jll|) and the KTeV and NA48 data 
we obtain 



Re e/e' = (1.9 ±0.46) x 10 



-3 



(2.4) 



where the error has been inflated according to the Particle Data Group procedure (a — > a x ^/% 2 /3) , to be used when 
averaging over exp erim ental data — in our case four sets — with substantially different central values. 

The value in eq. (2.4) can be considered the current experimental result. Such a result will be tested within the next 
year by the full data analysis from KTeV and NA48 and (hopefully) the first data from KLOE at DA$NE (Frascati); 
at that time, the experimental uncertainty will be reduced to a few parts in 10 4 . 

The most important outcome of the 1999 results is that direct CP violation has been unambiguously observed and 
that the superweak scenario, in which e' = 0, can be excluded to a high degree of confidence (more than 4 ct's). 
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III. HADRONIC MATRIX ELEMENTS 



In the present analysis, we use hadronic matrix elements for all the relevant operators Qi_io and a parameter Bk 
computed in the xQM at 0(p 4 ) in the chiral expansion jl7). This approach has three model-dependent parameters 
which are fixed by means of a fit of the AI — 1/2 rule. Let us here review briefly the model and how the matrix 
elements are computed. 

The xQM @1 is an effective quark model of QCD which can be derived in the framework of the extended Nambu- 
Jona-Lasinio model of chiral symmetry breaking (for a review, see, e.g., ]T9[]). 

In the xQM an effective interaction between the u, d, s quarks and the meson octet is introduced via the term 

£ X QM = —M (q R Y>q L +q L E T to) , (3.1) 

which is added to an effective low-energy QCD lagrangian whose dynamical degrees of freedom are the u, d, s quarks 



propagating in a soft gluon background. The matrix £ in (3.1) is the same as that used in chiral perturbation theory 



and it contains the pseudo-scalar meson multiplet. The quantity M is interpreted as the constituent quark mass in 
mesons (current quark masses are also included in the effective lagrangian). 

In the factorization approximation, the matrix elements of the four quark operators are written in terms of bet- 
ter known quantities like quark currents and densities. Such matrix elements (building blocks) like the current 
matrix elements (0| sj^ (1 — 75) u \K + (k)) and (7r + (p+)| s7 M (1 — 75) d \K + {k)) and the matrix elements of densities, 
(0| S75W \K + (k)), (tt + (p + )\ sd \K + (k)), are evaluated up to 0(p 4 ) within the model. The model dependence in the 
color singlet current and density matrix elements appears (via the M parameter) beyond the leading order in the 
momentum expansion, while the leading contributions agree with the well known expressions in terms of the meson 
decay constants and masses. 

Non-factorizable contributions due to soft gluonic corrections are included by using Fierz-transformations and by 
calculating building block matrix elements involving the color matrix T a : 

(0\s^T a (l- 75 )u\K+(k)) , (ir+(p+)\srT a (l- 75 )d\K + (k)) . (3.2) 

Such matrix elements are non-zero for emission of gluons. In contrast to the color singlet matrix elements above, they 
are model dependent starting with the leading order. Taking products of two such matrix elements and using the 
relation 

2 

9 2 s G%G a a p = \(^GG) (S^ - S^5 ua ) (3.3) 

makes it possible to express gluonic corrections in terms of the gluonic vacuum condensate pQ |. While the factoriz- 
able corrections are re-absorbed in the renormalization of the chiral couplings, non-factorizable contributions affect 
explicitly the form of the matrix elements. The model thus parameterizes all amplitudes in terms of the quantities 
M, (qq) , and (a s GG/n) . Higher order gluon condensates are omi tted . 



The hadronic matrix elements of the operatorial basis in eq. (1.5) for K — > kit decays, have been calculated 
up to 0{p A ) inclusive of chiral loops Jl?]]. The leading order (LO) (0(p°,p 2 )) matrix elements (Qi)f° and the 
NLO (0(p 2 ,p 4 )) corrections (Qi)f LO for final state isospin projections / = 0, 2 are obtained by properly combining 
factorizable and non-factorizable contributions and expanding the result at the given order. The total hadronic matrix 
elements up to 0(p 4 ) can then be written as: 

(Qi(M)>j = Z„y/Z^ [{Q^° + <Q t )f LO ] (ji) + al(p) , (3.4) 



where Qi are the operators in eq. (1.5), and aj(/j,) are the contributions from chiral loops (which include the wave- 
function renormalization). The scale dependence of the (Qi)j° :NLO comes from the perturbative running of the quark 
condensate and masses, the latter appearing explicitly in the NLO corrections. The quantities a\ (/i) represent the 
scale dependent meson- loop corrections which depend on the chiral quark model via the tree level chiral coefficients. 
They have been included to 0(p i ) in ref. |L7| by applying the MS scheme in dimensional regularization, as for the 
xQM calculation of the tree-level chiral coefficients. The wave-function renormalizations Zk and Z v arise in the xQM 
from direct calculation of the K — > K and 7r — > 7r propagators. 

The hadronic matrix elements are matched — by taking p-sd = fJ-LD — with the NLO Wilson coefficients at the scale 
A x ~ 0.8 (~ m p ) as the best compromise between the range of validity of chiral perturbation and that of strong 
coupling expansion. The scale dependence of the amplitudes is gauged by varying fi between 0.8 and 1 GeV. In this 
range the scale dependence of e'/e remains always below 10%, thus giving a stable prediction. 
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A. The fit of the AI = 1/2 rule 



In order to assign the values of the model-dependent parameters M, (qq) and (a s GG/'K) , we consider the CP- 
conserving amplitudes in the AI = 1/2 selection rule of K — ► tttt decays. In practice, we compute the amplitudes 



V2 cos d . 

G 1 



within the xQM approach and vary the parameters in order to reproduce their experimental values 
A (K -> tttt) = 3.3 x 10~ 7 GeV and A 2 (-K" -> tttt) = 1.5 x 1(T 8 GeV , 



(3.5) 
(3.6) 

(3.7) 



This procedure combines a model for low-energy QCD — which allow us to compute all hadronic matrix elements 
in terms of a few basic parameters — with the phenomenological determination of such parameters. In this way, some 
shortcomings of such a naive model (in particular, the matching between long- and short-distance components) are 
absorbed in the phenomenological fit. As a check, we eventually verify the stability of the computed observables 
against renormalization scale and scheme dependence. 

The fit of the CP-conserving involves the determination of the FSI phases. The absorptive component of the 
hadronic matrix elements appear when chiral loops are included. In our approach the direct determination of the 
rescattering phases gives at 0(p 4 ) So — 20° and 82 — —12°. Although these results show features which are in 
qualitative agreement with the phases extracted from pion-nucleon scattering [pl| , 



5 = 34.2° ±2.2° , 



8 2 = -6.9° ±0.2° 



(3.8) 



the deviation from the experimental data is sizeable, especially in the 1 = component. On the other hand, at 0(p A ) 
the absorptive parts of the amplitudes are determined only at 0(p 2 ) and disagreement with the measured phases 
should be expected. As a matter of fact, the authors of ref. [Q find that at 0(p 6 ) the absorptive part of the hadronic 
matrix elements are substantially modified to give values of the rescattering phases quite close to those in eq. (f3~S|)- 
At the same time the 0(p e ) corrections to the dispersive part of the hadronic matrix elements are very small. 

This result corroborates our ansatz |Q of trusting the dispersive parts of the 0(p A ) matrix elements while inputting 
the experimental values of the rescattering phases in all parts of our analysis, which amounts to taking cos<5o ss 0.8 
and cos 82 ~ 1. 
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1. Dependence of A$ and A2 on {qq}, (GG), Aq^ d 



and M at fi — 0.8 GeV. The gray and black sets of lines correspond 
to the extreme values of Aqcd and M. The length of the lines represents the effect of varying (qq), while keeping all other 
parameters fixed. The small dependence of A2 on the quark condensate is due to the the contribution of the electroweak 
penguins Q7,s- The vertical spread corresponds to varying (a s GG/n) , with the central line corresponding to the central value 
of (GG). The gray area denotes the region spanned by varying all the parameters without correlations in a ± 30% box around 
the experimental values of ^0 and A2 given by the cross hairs. The figure on the left (right) shows the HV (NDR) results, 
corresponding to varying {qq), {GG), Aqq D and M in the ranges given in eqs. ( |3.9| )-( ^ll[ ). 
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Hadronic matrix elements in the xQM depend on the 75-scheme utilized |L7j|. Their dependence partially cancels 
that of the short-distance NLO Wilson coefficients. Because this compensation is only numerical, and not analytical, 
we take it as part of our phenomenological approach. A formal 75-scheme matching can only come from a model more 
complete than the xQM. Nevertheless, the result, as shown in Fig. || below, is rather convincing. 

By taking 



A. 



(4) 
QCD 



340 ± 40 MeV 



(3.9) 



and fitting at the scale fi — 0.8 GeV the amplitudes in eqs. ( J3 . 5| ) — (3.6) to their experimental values, allowing for a ± 
30% systematic uncertainty, we find (see Fig. [l]) 



M = 195±?g MeV , 



(a s GG/ir) = (330 ± 5 MeV) 4 , {qq) = (-235 ± 25 MeV)' 



(3.10) 



in the HV-scheme, and 



M 



195t^ MeV : 



{a s GG/ir) 



(333±£ MeV) 



{qq) = (-245 ± 15 MeV)' 



(3.11) 



in the NDR-scheme. 

As shown by the light (NDR) and dark (HV) curves in Fig. 0, the 75-scheme dependence is controlled by the value of 
M, the range of which is fixed thereby. The 75 scheme dependence of both amplitudes is minimized for M ~ 190 — 200 
MeV. The good 75-scheme stability is also shown by e'/e and Bk- For this reason in our previous papers [[l7| we only 
quoted for these observables the HV results. 
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FIG. 2. Test of the 75-scheme stability of the Aq and A2 amplitudes as functions of M. The light (dark) curves corre- 
spond to the NDR (HV) results, while the horizontal lines mark the experimental values. T he tw o figures on the right plot 
AAi = 2\(A? v - A? DR )/{A? v + A? DR )\. For the values of (qq) and (a s GG/%) in eqs. (pTo|)-(^ll|) we find 75-scheme 
independence for M ~ 190 — 200 MeV. 



The fit of the amplitude Ao and A2 is obtained for values of the quark and gluon condensates which are in 
agreement with those found in other approaches, i.e. QCD sum rules and lattice, although it is fair to say that the 
relation between the gluon condensate of QCD sum rules and lattice and that of the xQM is far from obvious. The 
value of the constituent quark mass M is in good agreement with that found by fitting radiative kaon decays ]22| . 

In Fig. U we present the anatomy of the relevant operator contributions to the CP conserving amplitudes. It is 
worth noticing that, because of the NLO enhancement of the 1 = matrix elements (mainly due to the chiral loops), 
the gluon penguin contribution to Aq amounts to about 20% of the amplitude. 
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FIG. 3. Operator-by-operator contributions to Aq x 10 r and A2 x 10 8 (GeV) in the HV scheme. In (light) dark the (LO) 
NLO value. Notice the 0(p 4 ) enhancement of the gluon penguin operator in the I — amplitude. 



Turning now to the AS = 2 lagrangian, 



where 



T F" l, W 

4tt 2 



[\lr)xS(x c ) + \ 2 t -q 2 S{x t ) + 2X c X t r] 3 S(x c , x t )] 6(/j) 



where Xj — VjdV* s , Xi = rn^jm^. We denote by Qs2 the AS — 2 local four quark operator 

Qs2 = {sLl^d L )(s L l^d L ) , 



(3.12) 



(3.13) 



(3.14) 



which is the only local operator of dimension six in the standard model. 

The integration of the electroweak loops leads to the Inami-Lim functions (23) S(x) and S(x c , x t ), the exact expres- 
sions of which can be found in the reference quoted, which depend on the masses of the charm and top quarks and 
describe the AS = 2 transition amplitude in the absence of strong interactions. 

The short-distance QCD corrections are encoded in the coefficients 771, 772 and 773 with a common scale- and 
renormalization-scheme- dependent factor 6(/x) factorized out. They are functions of the heavy quarks masses and of 
the scale parameter Aqcd- These QCD corrections are available at the NLO |33| in the strong and electromagnetic 
couplings. 

The scale-dependent factor of the short-distance corrections is given by 



= [a. („)r 2/9 



(3.15) 



where J3 depends on the 75-scheme used in the regularization. The NDR and HV scheme yield, respectively: 



J. 



NDR 



3 162 3 162 

On the long-distance side, the hadronic Bk parameter is introduced by writing the AS* — 2 matrix element as 

4 



(K°\QS2^)\K°) = -f K mlB K (fi) 



(3.16) 



(3.17) 



The scale- and renormalization-scheme- independent parameter Bk is then defined by means of eqs. (3.15)— (3.17) as 

B K = b(jj)B K Qi). (3.18) 
By using the input values found by fitting the AI = 1/2 rule we obtain in both 75-schemes 

Br: = 1-1 ±0.2. (3.19) 



The result ( 3.19 ) includes chiral corrections up to 0(p 4 ) and it agrees with what we found in p7| . In the chiral limit 
one derives a simple expression for Bk (eq. (6.3) in ref. ]l7j]), which depends crucially on the value of the gluon 
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condensate. In this limit, for central values of the parameters, we obtain Bk = 0.44 (0.46) in the HV (NDR) scheme. 
A recent calculation of Bk, based on QCD sum rules, finds in the chiral limit and to the NLO in the 1/N expansion 
Bk = 0.41 ± 0.09 jH- Many calculations of Bk have been performed on the lattice (for a recent review see p5|). 
According to ref. |25f | , the analysis of ref. |^6]] presents the most extensive study of systematic errors and gives the 
(quenched) value Bk = 0.86 ± 0.06 ± 0.14, which should be taken as the present reference value, while awaiting for 
further progress in including dynamical quarks on the lattice. 

Notice that no estimate of e'/e can be considered complete unless it also gives a value for Bk- The case of the xQM, 
for instance, is telling insofar as the enhancement of Bq is partially compensated for by a large Bk (and accordingly 
a smaller Im At ) . 



B. The factors Bi 

The factors Bi, defined as 

B t = <Q s :) raodol /(Q,) VSA (3.20) 

have become a standard way of displaying the values of the hadronic matrix elements in order to compare them among 
various approaches. However they must be used with care because of their dependence on the renormalization scheme 
and scale, as well as on the choice of the VSA parameters. 




£9 

D 10 

B\ 2) ~ B< 2 > 

}(2) _ R (2) 



TABLE I. Central values of the Bi factors in the HV and NDR renormalization schemes. For B^^ the leading scaling 
dependence on (qq) is explicitly shown for a conventional value of the condensate. All othe r Bj fa ctors are either independent 
or very weakly dependent on (qq). The dependence on (a s GG/n) in the ranges of eqs. (3.10)-(3.11) remains always below 10%. 



They are given in the xQM m table [| in the HV and NDR schemes, at /i = 0.8 GeV, for the central value of Aq^ d . 
The dependence on Aqqd enters indirectly via the fit of the AI =1/2 selection rule and the determination of the 
parameters of the model. 

The uncertainty in the matrix elements of the penguin operators Q5-8 arises from the variation of (qq). This affects 
mostly the B^^ parameters because of the leading linear dependence on (qq) of the Q.5,6 matrix elements in the xQM, 
contrasted to the quadratic dependence of the corresponding VSA matrix elements. Accordingly, B 5 ^ scale as (qq) -1 , 
or via PCAC as m q , and therefore are sensitive to the value chosen for these parameters. For this reason, we have 
reported the corresponding values of B§ q when the quark condensate in the VSA is fixed to its PCAC value. 

It should however be stressed that such a dependence is not physical and is introduced by the arbitrary normalization 
on the VSA result. The estimate of e 1 is therefore almost independent of m q , which only enters the NLO corrections 
and the determination of Bk ■ 

The enhancement of the Q§fi matrix elements with respect to the VSA values (the conventional normalization of 
the VSA matrix elements corresponds to taking (qq)(0.8 GeV) ~ (—220 MeV) 3 ) is mainly due to the NLO chiral 
loop contributions. Such an enhancement, due to final state interactions, has been found in 1/A C analyses beyond 
LO p7|p8[ , as well as recent dispersive studies [EjlEjlJ. A large- N c approach, based on QCD sum rules |Q which 
reproduces the electroweak n + -Tr~ mass difference and the leptonic n (77) rare decays, disagrees in the determination 
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of {Qi)q 2 at /i = 0.8 GeV, due to the sharp scale dependence found for these matrix elements. Since the operator 
Q7 gives a negligible contribution to e'/e, we should wait for a calculation of other matrix elements within the same 
framework in order to asses the extent and the impact of the disagreement with the xQM results. 

Among the non-factorizable corrections, the gluon condensate contributions are most important for the CP- 
conserving 1 — 2 amplitudes (and account for the values and uncertainties of B\ 2) but are otherwise inessential 
in the determination of e', for which FSI are the most relevant corrections to the LO 1/N result. 



IV. BOUNDS ON Im At 



The updated measurements for the CKM elements iV^/Vc?,! implies a change in the determination of the Wolfenstein 
parameter 77 that enter in Im At. This is of particular relevance because it affects proportionally the value of e'/e. 
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FIG. 4. Bounds on the Wolfenstein parameters 77 = (1 — A 2 /2) and p = (1 — A 2 /2) and distribution of ImAt according to 
Parodi et al. for the input parameters in Table |n| 

The allowed values for ImAt — ^iKisllV^fcl 2 are found by imposing the experimental constraints for e, \V u b/V c b\, 
Anid and Am s which give rise to the following equations: 



where 



and 



A 2 
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The function F(x) in (4.5) in (11) can be found in Q. 



parameter 


value 




v ud 


0.9753 






0.2205 ± 0.0018 




sin t?iy 


0.2247 




mz 


91.187 GeV 






80.22 GeV 






4.4 GeV 




t ii ( rn 1 


1.3 GeV 


> 


\e\ 


(2.280 ±0.019) x 10 -3 


> 


\V cb \ 


0.0405 ± 0.0015 


> 


\v ub /v cb \ 


0.080 ±0.017 (CLEO), 0.104 ±0.019 (LEP) 




m t (mt) 


165 ± 5 GeV 




B K 


1.1 ±0.2 


> 


Am a 


> 12.4 ps" 1 


> 


Am d 


(0.472 ± 0.016) ps -1 


> 


J B d \J Bs d 


(2IO+32) MeV 


[> 


f — fr, /R„ 1 f„ / R„ 
? = JS S V DB B J jB d DB d 


1 1 1 +0.06 
r.ll_ .04 


> 




5.2792 ± 0.0018 GeV 


[> 




5.3693 ± 0.0020 GeV 


r> 




0.55 ±0.01 




f — f j. 


92.4 MeV 




f f 

J if — J K+ 


113 MeV 




'i t i — / i f i , f i-i n ] n 


138 MeV 




mx = m^o 


498 MeV 




m, 


548 MeV 




Ax 


2^2 TT.U 


> 




U.zo ± U.1U 






340 ± 40 MeV 




S. ± m d (1 GeV) 


12 ± 2.5 MeV 


> 


m 3 (1 GeV) 


150 ± 25 MeV 


> 




HV: (-235 ± 25 MeV) 3 , NDR: (-245 ± 15 MeV) 3 
HV: (330 ± 5 MeV) 4 , NDR: (3331^ MeV) 4 
HV: 195+^ MeV, NDR: 195+^ MeV 


> 
> 


(a s GG/7r) 
M 



TABLE II. Numerical values of the input parameters used in the present analysis. The triangles mark those that have been 
updated with respect to our 1997 estimate. 



By using the method of Parodi, Roudeau and Stocchi |34|j , who have run their program starting from the inputs 
listed in table [n] and 

r)i = 1.44 ± 0.18 , 772 = 0.52, r)a = 0.45 ± 0.01 , (4.7) 

which are the values we find for our inputs, it is found that 

ImA t = (1.14 ±0.11) x 10 -4 , (4.8) 

where the error is determined by the Gaussian distribution in Fig. [|. Notice that the value thus found is roughly 10% 
smaller than those found in other estimates for which Bk is smaller. 

The effect of this updated fit is a substantial reduction in the range of Im At with respect to what we used in our 
1997 estimate [l7|| : all values smaller than 1.0 x 10~ 4 are now excluded (as opposed as before when values as small as 
0.6 x 10~ 4 were included). 

V. ESTIMATING e'/e 

The value of e' computed by taking all input parameters at their central values (Table II) is shown in Fig. [| The 
figure shows the contribution to s' of the various operators in two 75 renormalization schemes at /1 = 0.8 GeV and 
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1.0 GeV. The advantage of such a histogram is that, contrary to the Bi, the size of the individual contributions does 
not depend on some conventional normalization. 
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Q3 Q4 Q5 Q6 Q7 Q8 Q9 Q10 All 
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FIG. 5. Contribution to e'/e (in units of 10 3 ) of each penguin operator in the HV and NDR 75-schemes. The figure on the 
left (right) corresponds to /1 = 0.8 (1.0) GeV. 



As the histogram makes it clear, the two dominant contributions come from the gluon and electroweak penguins, 
Qq and Qg. However, the gluon penguin dominates and there is very little cancellation with the electroweak penguin 
operator. The dominance of the 1 = components in the xQM originates from the 0(p 4 ) chiral corrections, the 
detailed size of which is determined by the fit of the AI = 1/2 rule. It is a nice feature of the approach that the 
renormalization scheme stability imposed on the CP conserving amplitudes is numerically preserved in e'/e. The 
comparison of the two figures shows also the remarkable renormalization scale stability of the central value once the 
perturbative running of the quark masses and the quark condensate is taken into account. 

In what follows, the model-dependent parameters M, (a s GG/ir) and (qq) are uniformly varied in their given ranges 
(flat scanning), while the others are sampled according to their normal distributions (see Table II for the ranges used). 
Values of e'/e found in the HV and NDR schemes are included with equal weight. 

0.8 1.3 1.7 2.1 2.5 2.9 3.3 3.4 3.8 4.2 4.6 



n 



Skewness = 0.9 
Median = 2.1 
Mean = 2.2, cr = 0.8 



n^n 
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0.04 
0.02 

y= 

0.8 1.3 1.7 2.1 2.5 2.9 3.3 3.4 3.8 4.2 4.6 

FIG. 6. Distribution of values of e'/e (in units of 10 -3 ). Normalized bins are plotted against the values of e'/e of each bin. 

For a given set, a distribution is obtained by collecting the values of e'/e in bins of a given range. This is shown in 
Fig. ^| for a particular choice of bins. The final distribution is partially skewed, with more values closer to the lower 
end but a longer tail toward larger values. However, because the skewness of the distribution is less than one, the 
mean and the standard deviation are a good estimate of the central value and the dispersion of values around it. This 
statistical analysis yields 

e'/e = (2.2 ±0.8) x 10~ 3 . (5.1) 

A more conservative estimate of the uncertainties is obtained via the flat scanning of the input parameters, which 
gives 

0.9 x 10~ 3 < Re e'/e < 4.8 x 10~ 3 . (5.2) 
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In both estimates a theoretical systematic error of ±30% is included in the fit of the CP conserving amplitudes A$ 
and A 2 . 

The stability of the numerical outcomes is only marginally affected by shifts in the value of fl v + V ' due to NLO chiral 
corrections |3q| and by additional isospin breaking effects |3^-|38| . Any effective variation of 0,,+,,' is anti-correlated 
to the value of (a s GG/ii) obtained in the fit of A 2 . We have verified that this affects the calculation of E>k and the 
consequent determination of Im A* in such a way to compensate numerically in e /e the change of Qn+ri' ■ Waiting for 
a confident assessment of the NLO isospin violating effects in the K — > hit amplitudes, we have used for £l v +ri' the 
'LO' value quoted in Table O. 

The weak dependence on some poorly-known parameters is a welcome outcome of the correlation among hadronic 
matrix elements enforced in our semi-phenomcnological approach by the fit of the AI = 1/2 rule. 

We have changed the central value of m c from 1.4 GeV, the value we used in |17|| , to 1.3 GeV in order to make our 
estimate more homogeneous with others. This change affects the determination of the ranges of the model parameters, 
mainly increasing (qq), via the fit of the CP conserving amplitudes. The central value of e'/e turns out to be affected 
below 10%. 

While the \QM approach to the hadronic matrix elements relevant in the computation of e'/e has many advantages 
over other techniques and has proved its value in the prediction of what has been then found in the experiments, it 
has a severe short-coming insofar as the matching scale has to be kept low, around 1 GeV and therefore the Wilson 
coefficients have to be run down at scales that are at the limit of the applicability of the renormalization-group 
equations. Moreover, the matching itself suffers of ambiguities that have not been completely solved. For these 
reasons we have insisted all along that the approach is semi-phenomenological and that the above shortcomings are 
to be absorbed in the values of the input parameters on which the fit to the CP conserving amplitudes is based. 



VI. OTHER ESTIMATES 



Figure [l] summarizes the present status of theory versus experiment. In addition to our improved calculation (and 
an independent estimate similarly based on the xQM), we have reported five estimates of e'/e published in the last 
year. Trieste's, Mimchen's and Roma's ranges are updates of their respective older estimates, while the other estimates 
are altogether new. 



X 2 
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Dortmund 



I 
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I 
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Lund 



Trieste 



FIG. 7. Theory vs. experiment in the year 2000. The gray/yellow band is the average experimental result. See the text for 
details on the various estimates. 



The estimates reported come from the following approaches: 

• Miinchcn's fl39| : In the Munchen approach (phenomenological l/N) some of the matrix elements are obtained 
by fitting the AI =1/2 rule at fj, = m c = 1.3 GeV. The relevant gluonic and electroweak penguin (Qq) and 
{Q%) 2 remain undetermined and are taken around their leading 1/iV values (which implies a scheme dependent 
result). In Fig. the HV (left) and NDR (right) results are shown. The darker range represents the result of 
Gaussian treatment of the input parameters compared to flat scanning (complete range) . 

• Roma's fjc|] : Lattice cannot provide us at present with reliable calculations of the 1 = penguin operators 
relevant to e'/e as well as of the 1 = components of the hadronic matrix elements of the current-current 
operators (penguin contractions), which are relevant to the AI =1/2 rule. This is due to large renormalization 
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uncertainties, partly related to the breaking of chiral symmetry on the lattice. In this respect, very promising 
is the domain-wall fermion approach |4lf| which allows us to decouple the chiral symmetry from the continuum 
limit. On the other hand, present lattice calculations compute K — > it matrix elements and use lowest order 
chiral perturbation theory to estimate the K — ► ttt: amplitude, which introduces additional (and potentially 
large) uncertainties. In the recent Roma re-evaluation of e' /e, B§ is taken at the VSA result varied by a 100% 
error. The estimate quotes the values obtained in two 75 schemes (HV and NDR). The dark (light) ranges 
correspond to Gaussian (flat) scan of the input parameters. 

Dortmund's p^j : In recent years the Dortmund group has revived and improved the approach of Bardeen, Buras 
and Gerard |42|| based on the 1/N expansion. Chiral loops are regularized via a cutoff and the amplitudes are 
arranged in a p 2n /N expansion. A particular attention has been given to the matching procedure between the 
scale dependence of the chiral loops and that arising from the short-distance analysis. The renormalization 
scheme dependence remains and it is included in the final uncertainty. The AI =1/2 rule is reproduced, but 
the presence of the quadratic cutoff induces a matching scale instability (which is very large for B^). The NLO 
corrections to (Qq) induce a substantial enhancement of the matrix element (right range in Fig. |7|) compared 
to the leading order result (left). The dark range is drawn for central values of m s , ^+17', IniA t and Aqcd- 



Dubna's 43 1: The hadronic matrix elements are computed in the ENJL framework including chiral loops up 
to Q(p 6 ) and the effects of scalar, vector and axial-vector resonances. (-Br-, and therefore ImA t , is taken 
from |n|) Chiral loops are regularized via the heat-kernel method, which leaves unsolved the problem of the 
renormalization scheme dependence. A phenomenological fit of the AI =1/2 rule implies deviations up to a 
factor two on the calculated (Qe). The reduced (dark) range in Fig. fj] corresponds to taking the central values 
of the NLO chiral couplings and varying the short-distance parameters. 

Taipei's J44| : Generalized factorization represents an attempt to parametrize the hadronic matrix elements in 
the framework of factorization without a-priori assumptions p5[ ] . Phenomenological parameters are introduced 
to account for non-factorizable effects. Experimental data are used in order to extract as much information as 
pos sible on the non-factorizable parameters. This approach has been applied to the K — > nn amplitudes in ref. 
[Q. The effective Wilson coefficients, which include the pcrturbative QCD running of the quark operators, are 
matched to the factorized matrix elements at the scale fj,p which is arbitrarily chosen in the perturbative regime. 
A residual scale dependence remains in the penguin matrix elements via the quark mass. The analysis shows 
that in order to reproduce the A J =1/2 rule and e'/e sizable non-factorizable contributions are required both 
in the current-current and the penguin matrix elements. However, some assumptions on the phenomenological 
parameters and ad hoc subtractions of scheme-dependent terms in the Wilson coefficients make the numerical 
results questionable. In addition, the quoted error does not include any short-distance uncertainty. 

Trieste's : The dark (light) ranges correspond to Gaussian (flat) scan of the input parameters. The bar on 
the left corresponds to the present estimate. That on the right is a new estimate }46| ], similarly based on the 
xQM hadronic matrix elements, in which however e'/e is estimated in a novel way by including the explicit 
computation of e in the ratio as opposed to the usual procedure of taking its value from the experiment. This 
approach has the advantage of being independent of the determination of the CKM parameters Im A* and of 
showing more directly the dependence on the long-distance parameter Bk as determined within the model. The 
difference (around 10%) between the two Tries te's estimates corresponds effectively to a larger value of ImAt, 



as determined from e only, with respect to eq. (4.8) 



Valencia's [£9|: The standard model estimate given by Pallante and Pich is obtained by applying the FSI 
correction factors obtained using a dispersive analysis a la Omnes-Mushkelishvili |^7| to the leading (factorized) 
1/N amplitudes. The detailed numerical outcome has been questioned on the basis of ambiguities related to the 
choice of the subtraction point at which the factorized amplitude is taken E] . Large corrections may also be 
induced by unknown local terms which are unaccounted for by the dispersive resummation of the leading chiral 
logs. Nevertheless, the analysis of ref. j2{| confirms the crucial role of higher order chiral corrections for e'/e, 
even though FSI effects alone leave the problem of reproducing the AI =1/2 selection rule open. 

Lund's [^9| : The AI = 1/2 rule and Bk have been studied in the NJL framework and 1/N expansion by 
Bijnens and Prades p8[ showing an impressive scale stability when including vector and axial- vector resonances. 
A recent calculation of e'/e at the NLO in 1/N has been performed in ref. |4S|| . The calculation is done in the 
chiral limit and it is eventually corrected by estimating the largest SU (3) breaking effects. Particular attention 
is devoted to the matching between long- and short-distance components by use of the A-boson method |30|,[5l| . 
The couplings of the A-bosons are computed within the ENJL model which improves the high-energy behavior. 
The AI = 1/2 rule is reproduced and the computed amplitudes show a satisfactory renormalization scale and 
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scheme stability. A sizeable enhancement of the Qg matrix element is found which brings the central value of 
e'/e at the level of 3 x 1(T 3 . 

Cut-off based approaches should also pay attention to higher-dimension operators which become relevant for match- 
ing scales below 2 GeV |p2| . The calculations based on dimensional regularization may avoid the problem if phe- 
nomenological input is used in order to encode in the hadronic matrix elements the physics at all scales. 

Other attempts to reproduce the measured e'/e using the linear a- model, which include the effect of a scalar 
resonance with m a ~ 900 MeV, obtain the needed enhancement of (Qq) ]5q |. However, the CP conserving 1 = 
amplitude falls short the experimental value by a factor of two. With a lighter scalar, m a ~ 600 MeV the CP 
conserving 1 = amplitude is reproduced, but e'/e turns out more than one order of magnitude larger than currently 
measured Q. 



VII. CONCLUSIONS 

The present analysis updates our 1997 estimate Re e'/e = 1.7 I^qX 10~ 3 |l7[], which already pointed out the potential 
relevance of non-factorizable contributions and the importance of addressing both CP conserving and violating data 
for a reliable estimate of e'/e. 

The increase in the central value is due to the update on the experimental inputs (mainly V u b/Vcb)- The uncertainty 
is reduced when using the Gaussian sampling, as opposed to the flat scan used in ref. |17[| . On the other hand the 
error obtained by flat scanning is larger due to the larger systematic uncertainty (±30%) used in the fit of the CP 
conserving amplitudes. 

Among the corrections to the leading 1/N (factorized) result, FSI play a crucial role in the determination of the 
gluon penguin matrix elements. Recent dispersive analyses of K — > rnr amplitudes show how a (partial) resummation 
of FSI increases substantially the size of the 1 = amplitudes, while slightly affecting the 1 = 2 components [p9[-pT{ . 
On the other hand, the precise size of the effect depends on boundary conditions of the factorized amplitudes which 
are not unambiguously known J4^,|55| . 

Finally, it is worth stressing that FSI by themselves do not account for the magnitude of the CP conserving decay 
amplitudes. In our approach a combination of non-factorizable soft-gluon effects (at 0(p 2 )) and FSI (at the NLO) 
makes possible to reproduce the Al = 1/2 selection rule. In turn, requiring the fit of the CP conserving K — > 7T7t 
decays allows for the determination of the "non-perturbativc" parameters of the xQM, which eventually leads to the 
detailed prediction of e'/e. Confidence in our final estimate of e'/e is based on the coherent picture of kaon physics 
which arises from the phenomenological determination of the model parameters and the self-contained calculation of 
all AS = 1 and 2 hadronic matrix elements. 
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